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1.  Introduction 


1.1  Background 

The  location  and  cleanup  of  buried  unexploded  ordnance  (UXO)  has  been  identified  as  a  high 
priority  mission-related  environmental  requirement  of  the  Department  of  Defense  (DoD).  The 
DoD  UXO  Response  Technology  Investment  Strategy  [1]  has  identified  wide  area  assessment  as 
one  of  six  technology  objectives,  with  a  goal  of  developing  capabilities  to  perform  rapid  initial 
assessment  of  large  areas.  The  benefit,  reiterated  by  the  recent  Defense  Science  Board  (DSB) 
Task  Force  on  UXO  [2]  is  that  the  identification  of  lands  where  UXO  is  not  present  may 
substantially  reduce  the  total  area  subjected  to  detailed  site  characterization  and  may  allow  for 
more  rapid  release  of  these  lands  to  reuse  functions. 

To  date,  a  great  deal  of  effort  in  the  UXO  research  community  has  been  focused  on  previously 
identified  sites  that  are  known  or  suspected  of  being  contaminated  with  UXO.  Total-coverage 
investigations  typically  involve  conducting  a  geophysical  survey  of  some  kind  (magnetometer, 
for  example)  with  sufficient  measurement  density  to  insure  that  measurements  are  taken  over  the 
entire  site.  The  results  are  then  analyzed  in  a  data-appropriate  manner  and  anomalies  are 
selected  and  prioritized.  Recommendations  for  remediation  are  made  based  on  the  results. 

Based  on  these  recommendations,  remediation  may  or  may  not  be  required  prior  to  land  reuse. 

The  broader  problem  posed  by  the  DSB  and  the  subject  of  the  ESTCP  Wide  Area  Assessment 
(WAA)  Pilot  Program  is  how  to  handle  large  areas  that  may  or  may  not  have  isolated  areas  of 
UXO  contamination  within  them.  Techniques  for  efficiently  surveying  a  large  area  to  locate 
concentrations  of  UXO  contamination  without  requiring  the  time-consuming  total-coverage  type 
of  survey  are  needed  to  address  this  problem.  We  have  demonstrated  one  such  technique  using  a 
ground-based  magnetometer  array  system. 

1.2  Objective  of  the  Demonstration 

A  comprehensive  WAA  program  has  the  possibility  of  making  an  immediate  impact  on  the  scope 
of  the  UXO  problem  currently  facing  the  DoD.  The  Defense  Science  Board  recently  estimated 
that  there  are  1400  sites  suspected  of  containing  UXO  contamination  covering  approximately  10 
million  acres  in  the  continental  US.  By  some  estimates,  as  much  as  80%  of  this  acreage  is  quite 
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likely  not  contaminated  with  UXO  at  all.  A  suite  of  technologies  that  can  accurately  and  rapidly 
delineate  the  areas  on  each  site  that  are  contaminated  from  those  that  are  not  contaminated  would 
lead  to  an  immediate  payback  in  terms  of  reducing  the  acreage  that  must  be  carefully  examined 
and  potentially  cleaned. 

We  have  demonstration  a  data  collection  and  analysis  methodology  to  support  the  rapid 
delineation  of  UXO  contamination  within  a  suspect  site.  Full-field  magnetometer  data  was 
collected  over  the  demonstration  site  along  planned  transects  provided  by  Pacific  Northwest 
National  Laboratories  (PNNL)  and  Sandia  National  Laboratory  (SNL)  in  cooperation  with  the 
ESTCP  Program  Office.  These  transects  were  designed  based  on  available  archive  information 
and  sound  statistical  sampling  methodologies.  Anomaly  location  and  a  measure  of  anomaly 
magnitude  were  extracted  from  these  data  using  an  automated  anomaly  detection  methodology. 
This  information  was  provided  to  PNNL  /  SNL  for  analysis  to  rapidly  delineate  UXO 
contamination  sites  such  as  impact  areas  and  bombing  targets.  With  the  rapid  pace  of  the 
automated  routines,  it  was  possible  to  interactively  plan  and  execute  additional  transects  to 
further  resolve  features  of  interest  while  the  survey  team  was  still  deployed  in  the  field. 


2.  Technology  Description 

2.1  Technology  Development  and  Application 

2.1.1  Vehicular  Magnetometer  System 

The  demonstration  was  conducted  using  the  Naval  Research  Laboratory  (NRL)  Multi-sensor 
Towed  Array  Detection  System  (MTADS).  The  MTADS  was  developed  with  support  from  the 
Environmental  Security  Technology  Certification  Program  (ESTCP).  The  MTADS  hardware 
consists  of  a  low-magnetic-signature  vehicle  that  is  used  to  tow  a  linear  array  of  eight 
magnetometer  sensors  over  large  areas  (25  acres  /  day)  to  detect  buried  UXO,  Eigure  1.  The 
sensors  are  sampled  at  50  Hz  and  typical  surveys  are  conducted  at  6  mph;  this  results  in  a 
sampling  density  of  ~6  cm  along  track  with  a  horizontal  sensor  spacing  of  25  cm.  Each 
magnetometer  measures  the  local  magnetic  field  of  the  earth  at  the  sensor. 

The  sensor  positions  are  measured  in  real-time  (5  Hz)  with  position  accuracies  of  '-'5  cm  using 
high  performance  Real  Time  Kinematic  (RTK)  Global  Positioning  System  (GPS)  receivers.  All 
navigation  and  sensor  data  are  time-stamped  with  Universal  Coordinated  Time  (UTC)  derived 
from  the  satellite  clocks  and  recorded  by  the  data  acquisition  computer  (DAQ)  in  the  tow 
vehicle.  The  positioning  technology  requires  the  availability  of  one  or  more  known  first-order 
survey  control  points.  The  sensor,  position,  and  timing  fdes  are  downloaded  periodically 
throughout  a  survey  onto  magnetic  disks  and  transferred  to  the  data  analyst  for  analysis. 

The  GPS  positioning  information  used  for  data  collection  is  shared  with  an  onboard  navigation 
guidance  display  and  provides  real-time  navigational  information  to  the  operator.  The  guidance 
display  was  originally  developed  for  the  airborne  adjunct  of  the  MTADS  system  (AMTADS)  [3] 
and  is  installed  in  the  vehicle  and  available  for  operator  use.  Eigure  2  shows  a  screenshot  of  the 
guidance  display  configured  for  vehicular  use.  An  integral  part  of  the  guidance  display  is  the 
ability  to  import  a  series  of  planned  survey  lines  (or  transects)  and  to  guide  the  operator  to  follow 
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these  transeets.  The  display  provides  left-right  course  correction  indicators,  an  optional  altitude 
indicator,  and  color-coded  flight  swath  overlays  where  the  current  transect  is  displayed  in  red 
and  the  other  transects  are  displayed  in  black  for  operator  reference.  The  survey  course-over¬ 
ground  (COG)  is  plotted  for  the  pilot  in  real  time  on  the  display  in  green.  If  a  GPS  outage  occurs 
during  the  survey,  the  COG  changes  color  to  warn  the  operator  and  allow  for  on-the-fly 
reacquisition  of  the  affected  area.  Figure  2  shows  the  operator  lining  up  to  survey  line  30  of  a 
transect  plan. 


Figure  1  -  MTADS  Magnetometer  system 


Figure  2  -  Screenshot  of  MTADS  Pilot  Guidance  Display 
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2,1,2  Data  Analysis  Methodology 

Each  data  set  is  collected  using  the  MagLogNT  software  package  (v2.921b,  Geometries,  Inc.). 
The  collected  raw  data  is  preprocessed  on  site  for  quality  assurance  purposes  using  standard 
MTADS  procedures  and  checks.  The  data  set  is  comprised  of  ten  separate  files,  each  containing 
the  data  from  a  single  system  device.  See  appendix  B  for  further  details  about  file  contents  and 
formats.  Each  device  has  a  unique  data  rate.  A  software  package  written  by  NRL  examines  each 
file  and  compares  the  number  entries  to  the  product  (total  survey  time  *  data  rate).  Any 
discrepancies  are  flagged  for  the  Data  Analyst  to  address.  Next,  the  data  is  merged  and  imported 
into  a  single  Oasis  montaj  (Geosoft,  Inc.)  database  using  custom  scripts  developed  from  the 
original  MTADS  DAS  routines  which  have  been  extensively  validated.  An  example  of  a 
working  screen  from  Oasis  montaj  is  shown  in  Eigure  3.  As  part  of  the  import  process  any  data 
corresponding  to  a  magnetometer  outage,  a  GPS  outage,  or  a  vehicle  stop  /  reverse,  is  defaulted 
or  marked  to  not  be  further  processed.  Defaulted  data  is  not  deleted  and  can  be  recovered  at  a 
later  time  if  so  desired.  Any  long  wavelength  features  such  as  the  diurnal  variation  of  the  earth’s 
magnetic  field  and  large  scale  geology  are  fdtered  from  the  data  (demedianed). 


Eigure  3  -  Working  screen  in  Oasis  montaj™  of  data  preprocessing  work  flow 


For  the  transect  surveys,  the  demedianed  magnetometer  data  are  converted  to  analytic  signal.  A 
built-in  feature  of  Oasis  montaj  is  used  to  extract  peaks  above  a  given  threshold  from  a  grid  like 
that  of  the  analytic  signal.  The  analytic  signal  is  used  because  anomaly  features  which  are 
dipolar  (having  both  positive  and  negative  components)  in  the  demedianed  magnetometer  data 
are  monopolar  in  the  analytic  signal.  The  detected  anomaly  locations  along  with  the  analytic 
signal  strength  at  the  peak  of  the  anomaly  were  provided  daily  to  the  ESTCP  Program  Office, 
PNNE,  and  SNL  for  the  previous  day’s  survey  results.  The  down-sampled  transect  COG  (6-10 


4 


m  spacing)  was  also  provided  at  the  request  of  PNNL  /  SNL.  The  data  analysis  work  flow  is 
shown  pietorially  in  Figure  4.  Additional  details  on  the  methodology  and  its  development  are 
available  in  Appendix  A. 

For  the  more  typieal  total  (100%)  eoverage  surveys,  the  loeated  demedianed  magnetometer  data 
were  imported  into  the  MTADS  Data  Analysis  System  (DAS)  software  for  individual  anomaly 
seleetion  and  analysis.  In  the  ease  of  isolated  munitions  in  the  far  field  (i.e.  farther  from  the 
sensors  than  their  eharaeteristie  dimension)  the  DAS  employs  resident  physies-based  models  to 
determine  target  size,  position,  and  depth.  A  spreadsheet  eontaining  details  of  the  anomaly 
loeation  and  fit  parameters  is  provided  along  with  the  loeations  of  anomalies  above  background 
which  are  identified  by  the  operator  but  for  which  the  dipole  model  do  not  give  a  reasonable  fit. 
The  loeated  demedianed  magnetometer  data  are  also  provided  for  arehival  purposes. 
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Figure  4  -  Automatic  anomaly  detection  scheme.  Example  data  is  from  the  MTADS  Test  Field 
at  Blossom  Point,  MD.  Magnetometer  data  is  shown  on  a  ±30  nT  vertieal  seale.  Analytie  signal 
data  is  shown  on  a  ±100  nT/m  vertieal  seale. 

2,2  Previous  Testing  of  the  Technology 

The  performanee  of  the  vehicular  MTADS  has  been  demonstrated  at  several  seeded  and  live 
ranges  sites  over  the  last  deeade  [4-9].  The  MTADS  has  demonstrated  probabilities  of  deteetion 
of  95  to  97%  and  loeation  aceuraeies  of  better  than  15  em  with  the  magnetometer  system  [7]. 

The  vehieular  MTADS  has  been  seleeted  to  serve  as  the  ground  truth  for  several  ESTCP- 
supported  demonstrations  of  potential  wide  area  survey  systems  [3,10,1 1]. 

As  an  example  of  the  performanee  of  the  MTADS,  the  results  from  the  survey  of  the  Target  SI  at 
Isleta  Pueblo,  NM  [11]  are  discussed  here  briefly.  Eor  the  Isleta  demonstration,  a  portion  of  the 
site  was  blind  seeded  by  the  ESTCP  Program  Offiee  with  a  variety  of  inert  munitions.  A  total 
eoverage  survey  was  condueted  over  the  site.  The  anomaly  list  generated  by  the  MTADS  was 
then  submitted  to  a  neutral  third  party  for  independent  evaluation.  The  results  were 
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representative  of  the  past  performanee  of  the  MTADS  system.  Analyzed  anomalies  were 
elassified  into  6  priority  eategories  where  1  is  likely  UXO,  3  is  unlikely  UXO,  4  is  unlikely  a 
elutter  item,  and  6  is  likely  a  elutter  item.  The  probability  of  deteetion,  Pd,  and  the  eumulative 
alarm  rate  were  determined  for  ineluding  eaeh  sueeessive  eategory  (from  1  to  6).  Pd  is  the 
fraetion  of  emplaeed  items  deteeted  and  the  false  alarm  rate  is  given  as  pieks  per  heetare  not 
eorresponding  to  an  emplaeed  item.  For  the  emplaeed  items  at  this  demonstration,  89%  of  the 
emplaeed  items  (Pd  =  0.89)  were  deteeted  and  plaeed  in  the  first  three  eategories  with  a  False 
Alarm  Rate  (FAR)  of  7  /  heetare.  The  loeational  performanee  metrics  were  mean  errors  of  -1 
and  4  cm  for  easting  and  northing,  respectively,  with  a  standard  deviation  of  12  and  13  cm  for 
the  same.  As  demonstrated  previously,  there  was  no  improvement  in  detection  by  widening  the 
detection  radius  from  1.0  to  1.5  m.  The  detection  radius  defines  how  large  an  error  in  reported 
position  can  still  be  considered  a  detection  of  the  emplaced  item. 


Several  hundred  detected  anomalies  were  selected  for  remediation  to  determine  the  performance 
of  the  systems  involved  in  the  overall  demonstration.  The  evaluation  metric  used  was  the 
location  difference  between  the  reported  location  of  the  anomaly  by  the  MTADS  and  the  actual 
location  reported  by  the  remediation  contractor.  As  was  seen  for  the  emplaced  items,  a  large 
majority  of  the  anomaly  picks  fall  well  within  the  more  restrictive  1.0-m  halo.  The  detailed 
location  performance  was  a  mean  miss  distance  of  35  cm.  90%  of  the  anomaly  picks  were 
within  59  cm  and  95%  were  within  77  cm  of  actual  remediated  location  of  the  anomaly.  As  was 
seen  for  the  emplaced  items,  a  large  fraction  of  the  remediated  anomalies  corresponding  to 
munitions  or  munitions-related  fragments  were  categorized  in  the  first  three  priority  groups  with 
95%  being  captured  in  the  first  two  priority  groups. 

2,3  Advantages  and  Limitations  of  the  Technology 

On  large  open  ranges  the  vehicular  MTADS  provides  an  efficient  survey  technology.  Surveys 
with  the  magnetometer  array  often  exceed  production  rates  of  20  acres  per  day.  UXO  items  with 
gauges  of  20mm  or  larger  are  typically  detected  to  their  self-burial  depths.  This  process  has  to 
date  involved  a  human  operator  manually  selected  the  data  corresponding  to  individual 
anomalies.  Each  data  segment  is  then  processed  by  a  physics-based  algorithm  incorporated  into 
the  MTADS  DAS  software. 

While  this  methodology  has  proven  highly  successful  in  the  past,  it  is  not  fast  enough  to  support 
the  rapid  data  requirements  for  the  transect  surveys  to  be  conducted  as  part  of  the  WAA  pilot 
project.  A  faster,  more  automated  method  was  developed  and  demonstrated.  The  location  and 
amplitude  of  anomalies  with  an  amplitude  above  an  empirically-determined  threshold  were 
reported  to  the  ESTCP  Program  Office,  PNNL,  and  SNL  along  with  the  survey  COG  for 
reference.  This  rapid  feedback  of  information  allowed  for  an  interactive  planning  and  execution 
of  additional  transects  while  the  demonstration  was  ongoing. 

The  presence  of  certain  terrain  features  such  as  deep  ravines  without  good  crossing  points  and 
thick  clusters  of  trees  limited  the  areas  that  could  be  surveyed  somewhat.  One  example  is  the 
“bowl”  feature  on  the  north  face  of  the  plateau  (south  of  Target  3,  Sections  3  and  4).  This  area 
was  so  broken  with  ravines  without  good  crossing  points  and  bound  by  the  plateau  on  the  south 
side  that  it  was  not  possible  to  reach  with  the  tow  vehicle.  In  other  cases,  the  presence  of  long 
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barbed-wire  fences  without  gates  and  deep  ravines  /  plateau  faces  without  good  access  points 
slowed  survey  operations  by  reducing  survey  line  length  and  increasing  travel  time  to  go  around 
these  obstacles. 


3.  Demonstration  Design 

3.1  Testing  and  Evaluation  Plan 

3.1.1  Demonstration  Set-Up  and  Start-Up 

The  former  Pueblo  Precision  Bombing  and  Pattern  Gunnery  Range  #2  (Pueblo  PBR#2)  is  located 
in  Otero  County,  Colorado,  approximately  20  miles  south  of  the  town  of  La  Junta  [12],  The 
training  range  encompasses  approximately  68,000  acres.  The  demonstration  area  encompasses 
approximately  7,400  acres  of  the  overall  Pueblo  PBR  #2  site  and  includes  Targets  3  and  4  along 
with  the  Suspected  75mm  Range  area  of  interest.  See  Reference  12  for  additional  discussion. 

The  coordinates  for  the  Pueblo  PBR#2  demonstration  site  are  given  in  Table  3-1. 


Table  3-1  -  Coordinates  for  the  Approximate  Comers  of  the  WAA  Pilot  Project  Pueblo  PBR  #2 
demonstration  site 


Point 

Latitude 

Longitude 

Northing  (m) 

Easting  (m) 

UTM  Zone  13N,  NAD  83 

SW 

37°  39'  52.662290656"  N 

103°  42'  02.32095666"  W 

4,169,400.00 

614,600.00 

MWl 

37°  42'  02.421430304"  N 

103°  42'  00.05663171"  W 

4,173,400.00 

614,600.00 

MW2 

37°  42' 01.969283698"  N 

103°  41'  19.22920282"  W 

4,173,400.00 

615,600.00 

NW 

37°  44'  44.16651 1803"  N 

103°  41'  16.36838703"  W 

4,178,400.00 

615,600.00 

NE 

37°  44'  42.784290086"  N 

103°  39'  13.81346694"  W 

4,178,400.00 

618,600.00 

MEl 

37°  43' 05.249205491"  N 

103°  39'  15.57919156"  W 

4,175,393.27 

618,600.00 

ME2 

37°  43' 05.360723355"  N 

103°  39'  25.35768396"  W 

4,175,393.27 

618,360.54 

SE 

37°  39'  50.892927079"  N 

103°  39'  24.40635276"  W 

4,169,400.00 

618,469.76 

The  MTADS  vehicular  system  mobilized  to  the  Pueblo  PBR#2  site  in  a  U.S.  Navy-owned  53-ft 
trailer.  The  tow  vehicle,  the  magnetometer  trailer,  notebook  computers  for  the  analysis  team, 
GPS  equipment,  batteries  and  chargers,  office  equipment,  radios  and  chargers,  tools,  equipment 
spares,  and  maintenance  items,  and  magnetometers  were  transported  in  the  trailer.  Harris 
Transportation  Company,  a  government  contract  transportation  firm  transported  the  trailer  to  the 
site. 

Due  to  the  remoteness  of  the  survey  site,  no  essential  support  services  were  available  on-site. 
Accordingly,  Nova  Research  made  provisions  to  acquire  all  of  the  requisite  supplies,  materials, 
and  facilities  from  local  rental  firms  to  establish  the  Pueblo  PBR#2  WAA  Base  Camp.  For  this 
operation  an  office  trailer  was  used  for  data  processing  and  analysis,  as  a  communications  center, 
for  battery  storage  and  charging  stations,  electronics  repair  station,  and  as  storage  for  spares  and 
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supplies.  This  trailer  was  supplied  with  AC  power,  heating,  and  eooling.  A  8’  x  40’  trailer, 
which  could  be  fully  opened  from  either  end  (for  drive-through  access),  was  provided  to  garage 
and  for  secure  storage  of  the  MTADS  vehicle  and  sensor  platform.  Power  to  the  trailers  was 
provided  by  a  diesel  field  generator  (65  kW  range)  that  was  also  used  to  recharge  the  vehicle, 
radios,  and  GPS  batteries  overnight.  Communications  among  on-site  personnel  was  provided  by 
hand-held  VHP  radios,  with  a  base  station  located  in  the  office  trailer.  Radios  were  provided  to 
all  field  and  office  personnel.  The  availability  of  cellular  phone  communications  on  site  was 
limited  but  was  available  at  the  higher  elevation  areas  of  the  site,  especially  in  the  general 
vicinity  of  GPS  control  points  Sky  CPS  and  Sky  CP6  and  at  the  U.S.  Forest  Service  corral  at  the 
comer  of  CR  25  and  CR  B.  Fuel  storage  was  provided  for  the  generator  with  refueling  as 
necessary  and  a  portable  toilet  was  provided  to  support  all  field  and  office  crews  with  weekly 
servicing.  Figure  5  shows  the  arrangement  of  this  logistics  support  at  a  recent  survey.  Due  to 
the  distance  from  the  WAA  Base  Camp  to  the  survey  areas  at  the  southern  end  of  the  site,  an 
additional  limited-scope  Auxiliary  Base  Camp  was  established  at  the  intersection  of  Roads  B  and 
23,  shown  in  Figure  6.  A  second  8’  x  40’  trailer  was  provided  to  garage  and  for  secure  storage  of 
the  MTADS  vehicle  and  sensor  platform  along  with  a  5  kW  generator  for  battery  charging. 

Upon  arrival  at  the  WAA  Base  Camp  on  August  28,  the  team  personnel  unpack  the  53’  trailer 
and  prepared  the  MTADS  vehicle  and  support  equipment  for  survey.  Another  performer  within 
the  WAA  Pilot  Project,  Sky  Research,  Inc.  has  established  eight  geodetic  survey  points  in  the 
general  area  of  the  demonstration  site.  The  coordinates  of  all  eight  points  are  given  in  Table  3-2 
(horizontal  datum:  North  American  Datum  of  1983,  1992  Adjustment  (NAD83/92);  vertical 
datum:  North  American  Vertical  Datum  of  1988  (NAVD88);  geoid  model:  National  Geodetic 
Survey  Geoid03). 


Figure  5  -  Photograph  of  the  WAA  Base  Camp  at  the  Pueblo  PBR#2  WAA  demonstration  site 
showing  the  relative  locations  of  the  trailers,  etc. 
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Figure  6  -  Photograph  of  the  Auxiliary  Base  Camp  at  Pueblo  PBR#2  WAA  demonstration  site 
Table  3-2  -  Survey  Control  Points  Installed  for  the  WAA  Pilot  Projeet  at  the  Pueblo  PBR  #2  site 


Point 

Name 

Latitude 

Longitude 

Ellipsoid 
Height  (m) 

Northing  (m) 

Easting  (m) 

Elevation  (m) 

WGS84 

UTM  Zone  13N,  NAD  83 

NAVD88 

Sky  CPI 

37°  39' 31.00828"  N 

103°  39'  23.98352"  W 

1397.643 

4168787.291 

618488.902 

1418.797 

Sky  CP2 

37°  39' 31.81861"  N 

103°  38' 50.55881"  W 

1393.904 

4168824.039 

619307.531 

1415.076 

Sky  CP3 

37°  38' 38.47452"  N 

103°  39' 56.81331"  W 

1401.753 

4167156.697 

617707.526 

1422.860 

Sky  CP4 

37°  39'  57.22970"  N 

103°  39'  23.93856"  W 

1396.392 

4169595.462 

618478.424 

1417.560 

Sky  CPS 

37°  44' 38.76102"  N 

103°  40'  42.21410"  W 

1439.133 

4178245.156 

616438.238 

1460.417 

Sky  CP6 

37°  44'  16.35566"  N 

103°  39'  18.00019"  W 

1453.184 

4177583.970 

618509.241 

1474.494 

Sky  CP7 

37°  42'  07.59478"  N 

103°  38'  13.14673"W 

1365.804 

4173638.425 

620154.391 

1387.077 

Sky  CPS 

37°  43' 01.02076"  N 

103°  42' 11.45205"  W 

1468.740 

4175202.188 

614295.936 

1489.925 

On  August  29,  the  RTK  GPS  base  station  receiver  and  radio  link  were  established  on  one  of  the 
available  established  control  points  (initially  Sky  CP6,  Sky  CP5  and  Sky  CPI  were  used  as 
required  during  the  demonstration).  The  magnetometer  trailer  was  connected  to  the  tow  vehicle 
and  the  system  was  powered  up.  The  connectivity  of  the  magnetometers  to  the  DAQ  computer 
and  the  establishment  of  normal  SNR  performance  were  verified  along  with  the  operational  state 
of  the  vehicle  RTK  GPS  system.  A  period  (5-6  minutes)  of  quiet,  static  data  was  collected  and 
submitted  to  the  Data  Analyst  for  validation.  This  same  static  test  was  repeated  throughout  the 
survey  campaign  each  morning  prior  to  surveying  the  calibration  items. 

A  lane  of  emplaced  calibrations  items  was  emplaced  by  the  demonstration  team  and 
representatives  of  the  ESTCP  Program  Office  south  of  the  WAA  Base  Camp  between  the  WAA 
Base  Camp  and  the  WAA  demonstration  area.  The  schedule  of  ground  based  system  calibration 
items  emplaced  for  the  Pueblo  PBR#2  site  is  given  in  Table  3-3.  A  multi-pass  magnetometer 
survey  of  the  proposed  calibration  strip  was  conducted  prior  to  emplacement  and  the  quietest 
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area  in  terms  of  geology  was  seleeted  for  the  ealibration  items.  The  eomposition  of  the  ground  in 
the  seleeted  area  was  approximately  0.5  m  of  soil  and  broken  roek  on  top  of  a  hard  rook  layer. 
Consequently,  emplaoement  depths  were  limited  to  60  om.  Once  the  items  were  emplaced  and 
photographed,  the  positions  of  each  item’s  nose  and  tail  were  recording  using  RTK  GPS.  The 
digital  photographs  of  the  emplaced  items  are  available  on  the  attached  DVD.  The  holes  were 
refilled  with  the  removed  material  and  leveled.  A  single  pass  magnetometer  survey  was 
conducted  over  the  emplaced  items.  Additional  surveys  of  the  calibration  items  were  conducted 
at  the  beginning  and  end  of  each  work  day.  On  a  few  occasions  it  was  not  possible  to  collect  the 
end-of-day  calibration  data  due  to  range  closure  due  to  weather  or  equipment  malfunction 
bringing  about  an  abrupt  end  to  the  day’s  efforts.  Prior  to  operating  out  of  the  Auxiliary  Base 
Camp,  Sphere  #1  was  relocated  to  near  the  Auxiliary  Base  Camp  for  calibration  purposes.  The 
new  location  of  the  sphere  is  also  given  in  Table  3-3. 

Table  3-3  -  Schedule  of  Ground-based  System  WAA  Calibration  Items  for  Pueblo  PBR#2 


ID 

UTM  Zone  13N 

Actual 

Easting  (m) 

Northing  (m) 

Depth  (cm) 

Grid  Azimuth 
(deg) 

North  Cali 

bration  Lane 

Sphere  #1  (Driver  Side) 

616434.500 

4178732.403 

0 

N/A 

Sphere  #2  (Passenger  Side) 

616435.459 

4178732.109 

0 

N/A 

155mm  Projectile  #2 

616441.180 

4178749.703 

35 

35 

60  mm  Mortar  #2 

616447.267 

4178768.742 

30 

46 

105mm  Projectile  #2 

616453.828 

4178787.621 

60 

44 

105mm  Projectile  #1 

616459.639 

4178806.967 

45 

178 

81mm  Mortar  #2 

616465.341 

4178825.771 

43 

69 

81mm  Mortar  #1 

616469.792 

4178839.941 

25 

20 

155mm  Projectile  #1 

616474.350 

4178854.198 

50 

46 

60  mm  Mortar  #1 

616478.260 

4178868.186 

10 

148 

37mm  Sim  #2 

616481.220 

4178877.825 

10 

57 

37mm  Sim  #1 

616484.096 

4178887.419 

5 

160 

South  Calibration  Sphere 

Sphere  #1  (Driver  Side) 

618349.507 

4168713.105 

0 

N/A 

The  Site  Safety  Officer  conducted  a  ‘tail-gate’  safety  meeting  prior  to  beginning  the  day’s  efforts 
and  again  each  day  that  personnel  were  on  site.  The  topic(s)  for  each  day’s  meeting  were  at  the 
discretion  of  the  Site  Safety  Officer  and  focused  on  safety  issues  relating  to  the  day’s  planned 
work. 

Preventative  maintenance  inspections  were  conducted  at  least  once  a  day  by  all  team  members, 
focusing  particularly  on  the  tow  vehicle  and  sensor  trailer.  Any  deficiencies  were  addressed 
according  to  the  severity  of  the  deficiency.  Parts,  tools,  and  materials  for  many  maintenance 
scenarios  are  available  in  the  system  spares  inventory  located  on  site.  Routine  tools  and 
supplies,  for  example  spare  tires  for  the  tow  vehicle  and  sensor  trailer,  were  carried  in  the  chase 
vehicle  which  accompanied  the  tow  vehicle  onto  the  site.  Status  on  break-downs  /  failures  that 
resulted  in  long-term  delays  in  surveying  was  reported  to  the  WAA  Project  Manager  as 
appropriate. 


10 


3.1,2  Period  of  Operation 


The  main  portion  of  the  demonstration  was  aecomplished  from  Tuesday,  August  30'*'  through 
Saturday,  Oetober  22"‘*.  Operations  were  eondueted  in  three  portions  as  detailed  in  tabular  form 
in  Table  3-4.  The  originally  scheduled  second  survey  portion  was  broken  into  two  due  to 
unscheduled  maintenance  required  on  the  tow  vehicle. 

Table  3-4  -  Pueblo  PBR  #2  Survey  Demonstration  Deployment  Schedule 


Date 

Planned  Action 

Week  of  August  15*'' 

Pack  trailer  at  Blossom  Point. 

Friday,  August  19"' 

Trailer  leaves  Blossom  Point  for  Pueblo  PBR  #2. 

Tuesday,  August  23**' 

Trailer  arrives  Pueblo  PBR  #2. 

Sun,  August  28*'' 

Personnel  arrive  La  Junta;  unpack  trailer,  assemble  MTADS 
system. 

Mon,  Aug  29**' 

Scout  demonstration  area,  emplace  and  survey  calibration  items. 

Tue,  Aug  30**' 

Begin  ground  surveys. 

Fri,  Sep  16**' 

Pause  ground  surveys.  Personnel  departs  site. 

Sun,  Oct  2**' 

Personnel  return  to  La  Junta. 

Mon,  Oct  3'" 

Resume  ground  surveys. 

Fri,  Oct  7**' 

Pause  ground  surveys,  arrange  for  vehicle  maintenance. 

Sat,  Oct  8**'  -  Tue,  Oct 

11**' 

Personnel  depart  La  Junta. 

Mon,  Oct  17**' 

Personnel  return  to  La  Junta,  reassemble  vehicle. 

Tue,  Oct  18**' 

Resume  ground  surveys. 

Sat,  Oct  22**' 

Complete  ground  surveys. 

Sun,  Oct  23**' 

Pack  trailer. 

Mon,  Oct  24**' 

Personnel  depart  La  Junta. 

Thu,  Nov  10**' 

Trailer  departed  Pueblo  PBR  #2. 

Mon,  Nov  14**' 

Trailer  arrives  at  Blossom  Point,  MD. 

Week  of  Nov  28**' 

Submit  Data  Report  to  ESTCP. 
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3.1,3  Operational  Parameters  for  the  Technology 

The  precision  collection  of  high  SNR  magnetometer  data  using  the  MTADS  platform  is  a  mature 
technology.  The  rapid  and  accurate  extraction  of  anomaly  location  and  a  measure  of  anomaly 
amplitude  (peak  analytic  signal)  from  high-volume  transect  data  collection  is  the  novel 
component  of  this  demonstration.  To  accomplish  this  task  an  automated  method  of  extracting 
the  anomaly  locations  from  the  survey  data  was  required.  One  such  method  has  been  developed 
and  is  discussed  in  detail  in  Appendix  A.  Briefly,  the  located  magnetic  field  data  (nT)  are 
collected  as  normal  for  an  MTADS  survey.  The  demedianed  total  field  data  are  converted  to 
analytic  signal  (AS,  nT/m)  where  the  analytic  signal  is  calculated  from  the  squares  of  the 
derivatives  in  the  x,  y,  and  z  directions: 


This  process  involves  a  gridding  step,  where  real-world  data  is  interpolated  onto  a  line-scale 
mesh  with  a  defined  grid  cell  size.  The  use  of  a  regular  grid  reduces  the  complexity  of  the 
calculations  required  for  the  following  steps.  The  utility  of  the  analytic  signal  is  that  anomaly 
features  which  are  dipolar  (have  both  positive  and  negative  components)  in  the  total  field  are 
monopolar  in  the  analytic  signal.  This  facilitates  the  detection  of  anomalies. 

One  can  then  define  the  peak  cut-off  threshold,  grid  smoothing  parameters  required  to  eliminate 
multiple  picks  per  target,  and  grid  cell  size  to  be  used  for  the  analysis.  Initial  analysis  (See 
Appendix  A)  has  shown  that  these  parameters  may  be  similar  for  several  sites  with  diverse 
geology  and  have  the  potential  to  be  applied  more  generally.  This  assertion  was  evaluated 
during  the  early  data  collection  stages  by  optimizing  two  parameters  (#  of  smoothing  passes,  and 
peak  threshold  cut-off  value)  against  the  incoming  data.  The  grid  cell  size  parameter  was  not 
varied  as  initial  testing  indicated  that  processing  times  become  prohibitive  at  grid  cell  sizes 
smaller  than  0.25m  for  transects  of  any  length.  When  the  survey  results  from  the  calibration  strip 
and  several  transect  data  sets  from  the  first  day  of  data  collection  were  available,  the  data  was 
used  to  evaluate  the  extraction  parameters  for  the  Pueblo  PBR#2  site.  At  a  fixed  cut-off 
threshold  of  25  nT/m,  the  effect  of  the  number  of  passes  of  the  smoothing  filter  was  evaluated 
and  the  initial  result  of  6  passes  was  found  appropriate.  Such  an  evaluation  on  an  early  data  set 
(05242006)  from  August  30**'  is  shown  in  Figure  7.  The  RMS  variation  in  the  analytic  signal 
from  quiet  portions  of  the  data  was  evaluated  and  found  to  be  on  the  order  of  2.5  nT/m.  Starting 
with  this  level  and  increasing  the  cut-off  threshold  in  increments  of  2.5  nT/m,  the  anomaly 
extraction  results  were  determined  and  a  peak  anomaly  cut-off  value  of  20-25  nT/  was  found  to 
be  necessary  to  avoid  extracting  spurious  anomalies  from  the  incoming  data.  The  results  for  an 
early  data  set  (05242006)  from  August  30**^  are  shown  in  Figure  8.  A  threshold  of  25  nT/m  was 
used  for  the  remainder  of  the  demonstration  transect  results. 
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Figure  7  -  Effect  of  increasing  number  of  passes  of  smoothing  filter  on  the  05242006  data  set 
results.  The  red  line  indicates  the  result  for  the  final  parameter  value. 


Figure  8  -  Effect  of  increasing  peak  anomaly  cut-off  threshold  value  on  the  05242006  data  set 
results.  The  red  line  indicates  the  result  for  the  final  parameter  value. 

3,1,4  Transect  Magnetometer  Survey  Results 

The  major  focus  of  the  data  collection  effort  for  the  first  three  weeks  of  the  demonstration  was 
the  collection  of  transect  magnetometer  data  following  the  transect  plans  provided  by  PNNL  / 
SNL  based  on  archive  data  (CSM  vO)  and  WAA  Pilot  Project  goals.  The  transect  plans  were 
divided  into  three  categories:  1)  Transects  to  interrogate  the  entire  WAA  Demonstration  site  for 
the  actual  positions  and  footprints  of  Targets  3  and  4  as  noted  in  CSM  vO  and  any  additional 
similar  features  of  interest,  2)  East  /  West  Transects  to  interrogate  the  Suspected  75mm  Range 
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area  of  interest  for  possible  features  of  interest,  and  3)  Any  additional  transeets  requested  by 
PNNL  /  SNL  /  ESTCP  Program  Offiee  based  on  the  results  of  items  1  &  2  or  other  data. 

For  the  first  eategory,  N/S  transeets  eovering  the  entire  WAA  demonstration  area,  two  transeet 
designs  were  prepared  by  PNNL/SNL.  The  first,  sparse  design  was  based  on  traversing  100-lb 
praetiee  bomb  targets  and  features  of  interest  with  a  99%  probability  of  traversing  a  1000  ft 
cireular  target  or  feature  of  interest.  The  transeets  were  oriented  N/S  with  a  308  m  spaeing.  The 
seeond,  eonservative  design  was  based  on  finding  500-ft  diameter,  cireular  100-lb  practice  bomb 
targets  with  a  99%  probability  of  traversing  the  target  or  feature  of  interest.  The  transects  were 
oriented  N/S  with  a  154  m  spacing.  This  design  leverages  the  data  already  collected  as  part  of 
the  sparse  design  and  adds  an  addition  transect  equally  spaced  between  each  pair  of  sparse 
transects. 

For  the  second  category,  F/W  transects  covering  the  suspected  75mm  range,  two  designs  were 
prepared  by  PNNF/SNF.  The  first,  sparse  design  was  based  on  a  99%  probability  of  traversing  a 
100  m  X  400  m  elliptical  target  or  feature  of  interest.  The  transects  were  oriented  F/W  with  400 
m  spacing  and  leveraging  the  N/S  transects  already  recorded.  The  second,  conservative  design 
was  based  on  a  99%  probability  of  traversing  a  100  m  diameter,  circular  feature  of  interest.  The 
transects  were  oriented  E/W  with  100  m  spacing  and  leveraging  the  N/S  transects  already 
recorded.  This  design  leverages  the  data  already  collected  as  part  of  the  sparse  design  and  adds 
three  additional  transect  equally  spaced  between  each  pair  of  sparse  E/W  transects. 

For  the  third  category,  17  100-m  spacing  E/W  transects  starting  from  50m  north  of  the  southern 
boundary  of  the  WAA  demonstration  site  were  surveyed  to  further  define  the  footprint  of  Target 
4.  Four  additional  areas  of  interest  were  also  identified  from  the  N/S  transect  data  by  PNNF  / 
SNF,  labeled  Areas  23,  25,  26,  27.  Based  on  CSM  vl,  four  additional  areas  of  interest  were 
identified.  In  these  additional  areas  of  interest,  transect  plans  of  4  -  10  transects  were  designed 
and  surveyed.  Oasis  montaj  .XYZ  files  of  all  planned  transects  are  included  on  the  attached 
DVD.  As  an  example,  a  portion  of  the  N/S  transect  plan  shown  in  Figure  9  along  with  the  COG 
of  the  transect  data  collected  on  September  5,  2005. 

The  position  (easting,  northing)  and  signal  strength  (peak  analytic  signal)  were  extracted  for  each 
anomaly  above  an  empirically  determined  threshold  for  all  transect  data.  Data  collection  began 
with  transects  that  did  not  included  identified  targets  to  allow  for  the  establishment  of  the  cut-off 
threshold  based  on  the  determined  noise  floor.  Figure  10  shows  the  results  of  all  transect  data 
collected  in  course  of  this  demonstration.  The  COGs  are  shown  as  green  lines  and  each  detected 
anomaly  is  shown  as  a  filled  circle. 

The  total  acreage  covered  by  transect  surveys  was  143  acres,  or  approximately  2%  of  the  total 
7,400  acres  site.  Natural  topology  (ravines,  plateau  faces,  trees,  etc.)  and  man-made  obstructions 
(e.g.  fences)  made  it  difficult  and  impractical  to  complete  each  transect  in  a  single  survey. 
Therefore  each  transect  was  broken  into  one  or  more  segments  in  the  field.  The  flexibility  of  the 
MTADS  Pilot  Guidance  software  allows  for  this  to  be  done  easily  and  on  the  fly.  The  exact 
details  of  the  area  covered  by  each  survey  file  are  given  in  an  Excel  spreadsheet  on  the  attached 
DVD  (Transects_Summary.xls).  An  excerpt  of  the  annotated  listing  is  given  in  Table  3-5.  The 
corresponding  demedianed  magnetometer  data,  the  analytic  signal  data,  the  anomaly  list,  and  the 
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COG  files  for  each  transect  survey  are  also  supplied  on  the  attached  DVD  in  the  “Transect 
Surveys”  subdirectory. 


a)  b) 

Figure  9  -  Map  showing  a)  North  /  South  Transect  Plan  (sparse  shown  in  red,  additional 
transects  for  conservative  shown  in  green),  b)  Transect  plan  with  actual  survey  for  Julian  date 
(05248,  September  5)  shown. 
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Figure  10  -  Map  showing  the  transect  survey  results  for  the  Pueblo  PBR#2  demonstration. 
Transect  COGs  are  shown  as  green  lines  and  individual  detected  anomalies  are  shown  as  filled 
circles. 
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Table  3-5  -  Excerpt  of  Annotated  Listing  of  Transect  Surveys  Conducted  During  Pueblo  PBR#2 
Demonstration. 


Day  1  Deliverables 

05242004 

Northern-most  portion  of  NS  Line  21 

05242006 

Northern-most  portion  of  NS  Line  19 

05242007 

Northern-most  portion  of  NS  Line  17 

05242008 

Northern-most  portion  of  NS  Line  15 

05242010 

Northern-most  portion  of  NS  Line  1 1 

05242011 

Northern-most  portion  of  NS  Line  9 

05242012 

Northern-most  portion  of  NS  Line  7 

05242013 

Northern-most  portion  of  NS  Line  12 

05242014 

Northern-most  portion  of  NS  Line  14 

05242015 

Northern-most  portion  of  NS  Line  16 

05242016 

Northern-most  portion  of  NS  Line  18 

Day  2  Deliverables 

05243003 

Northern-most  portion  of  NS  Line  8 

05243004 

Northern-most  portion  of  NS  Line  10 

05243005 

Northern-most  portion  of  NS  Line  13 

05243006 

Northern-most  portion  of  NS  Line  20 

05243007 

Northern-most  portion  of  NS  Line  22 

05243008 

Northern-most  portion  of  NS  Line  24 

05243009 

Northern-most  portion  of  NS  Line  26 

Day  3  Deliverables 

05244004 

Northern-most  portion  of  NS  Line  23 

05244005 

Northern-most  portion  of  NS  Line  27 

05244006 

NS  Line  7  (Plateau  in  Seetion  8) 

05244007 

NS  Line  8  (Plateau  in  Seetion  8) 

05244008 

NS  Line  9  (Plateau  in  Seetion  8) 

05244009 

NS  Line  27  (South  of  ravine  in  Seetion  2) 

05244010 

NS  Line  26  (South  of  ravine  in  Seetion  2) 

05244011 

NS  Line  24  (South  of  ravine  in  Seetion  3) 

05244012 

NS  Line  23  (South  of  ravine  in  Seetion  3) 

05244013 

NS  Line  22  (South  of  ravine  in  Seetion  3) 

05244014 

NS  Line  21  (South  of  ravine  in  Seetion  3) 

05244015 

NS  Line  20  (South  of  ravine  in  Seetion  3) 

05244017 

NS  Line  23  (Simmons  half-seetion) 

05244018 

NS  Line  22  (Simmons  half-seetion) 

05244019 

NS  Line  21  (Simmons  half-seetion) 

... 

... 

3.1.5  Total  Coverage  Magnetometer  Survey  Results 

In  addition  to  the  transect  surveys  covering  the  breadth  of  the  WAA  demonstration  area,  several 
small  areas  (30  -  85  acres)  were  selected  for  total  coverage  surveys.  Areas  were  selected  in 
cooperation  with  the  ESTCP  Program  Office  to  achieve  three  objectives:  1)  Collect  data  in  areas 
identified  by  the  transect  surveys  as  “quiet”  to  determine  what  the  background  anomaly  density 
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for  the  WAA  demonstration  site,  2)  Colleet  data  near  Targets  3  and  4  to  evaluate  the  anomaly 
density.  By  starting  near  the  target  and  moving  away  in  several  steps,  it  is  possible  to  map  the 
anomaly  density  falloff  as  one  moves  away  from  the  Target,  and  3)  Colleet  addition  data  on  the 
Suspeeted  75mm  Range  area  of  interest  in  support  of  the  transect  survey  results.  These  surveys 
were  conducted  as  typical  MTADS  magnetometer  surveys  with  a  line  spacing  of  2.0  m  (tire  next 
to  tire  spacing).  Collected  and  processed  magnetometer  data  were  exported  from  the  Oasis 
montaj  environment  and  loaded  into  the  MTADS  DAS  software  for  individual  anomaly  analysis. 
Detailed  target  lists  for  each  area  are  provided  on  the  attached  DVD  along  with  the 
magnetometer  data  archives.  See  Appendix  B  for  the  file  format  details. 

Figure  1 1  shows  the  total  coverage  area  anomaly  maps  superimposed  on  the  WAA 
demonstration  site  topographical  map.  Table  3-6  contains  a  summary  of  the  total  coverage 
survey  results.  Column  three  lists  the  number  of  anomalies  extracted  by  the  operator  in  the  DAS 
in  each  area  and  column  five  lists  the  number  of  those  anomalies  which  could  be  fit  using  the 
resident  dipole  model  to  a  coherence  value  of  0.85  or  better. 

Table  3-6  -  Total  Coverage  Area  Result  Summary 


Target 

Area 

Number  of 
Anomalies 

Anomalies  / 
Acres 

Number  of 
Dipole  Fits 

Acres 

Target  4 

BT4  Center 

938 

85 

873 

11.0 

1C 

1095 

28 

938 

38.8 

IB 

245 

7 

242 

33.9 

lA 

169 

5 

168 

33.6 

Target  3 

3A 

2112 

60 

1830 

35.4 

3B 

520 

14 

519 

36.3 

3C 

207 

6 

206 

35.7 

Simmons  Area 

72 

1 

72 

85.0 

Suspected  75mm 
Range 

2A 

148 

5 

148 

31.0 

2B 

83 

2 

83 

36.8 

An  85  acre  area  was  selected  in  the  northern  portion  of  Section  10,  referred  to  as  the  Simmons 
Area  in  reference  to  the  ranchers  who  currently  lease  the  area  (Brian  and  Janet  Simmons),  as  a 
“quiet”  area  for  determining  the  background  anomaly  level.  The  transect  survey  results  indicated 
that  this  area  had  very  few  anomalies  with  only  a  single  anomaly  detected  by  the  transect  surveys 
within  the  Simmons  Area  total  coverage  area.  Figure  12  presents  the  magnetometer  data 
anomaly  map  for  the  Simmons  Area. 

Four  total  coverage  areas  were  surveyed  in  the  vicinity  of  Target  4,  located  in  the  southern  part 
of  the  WAA  demonstration  site.  The  three  main  total  coverage  areas  are  labeled  Area  lA,  IB, 
and  1C  with  Area  1C  being  the  closest  to  Target  4  and  Area  lA  being  the  furthest  east.  The  area 
BT4  Center  was  part  of  an  earlier  survey  scheme  developed  to  map  the  anomaly  density  falloff 
from  Target  4  which  was  altered  to  the  Area  lA  -  1C  plan  at  the  request  of  the  landowners 
involved  (Ralph  and  Russell  Rounds).  The  BT4  Center  data  consists  of  4  acres  and  is  presented 
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for  completeness.  Figure  13,  Figure  14,  Figure  15,  and  Figure  16  present  the  magnetometer  data 
anomaly  maps  for  Areas  1C,  IB,  lA,  and  BT4  Center  respeetively. 

The  data  shown  in  the  magnetometer  anomaly  maps  shown  for  the  total  eoverage  areas  have 
been  masked  by  the  final  perimeters  of  the  areas.  The  are  hive  data  sets  are  also  bounded  by 
these  perimeters.  Oasis  montaj  polygon  files  (.ply)  of  eaeh  perimeter  are  provided  on  the 
attaehed  DVD.  Data  outside  the  perimeters  (turn  arounds,  ete.)  is  generally  unreliable  for  dipole 
analysis  beeause  of  severely  non-parallel  direetions  but  alternative  data  sets  with  all  data  are  also 
provided  on  the  attaehed  DVD.  These  data  should  be  used  with  eaution. 
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Figure  1 1  -  Pueblo  PBR#2  Total  Coverage  Survey  Results 
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Figure  12  -  Magnetometer  Anomaly  Map  of  the  Simmons  Area 
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Figure  13  -  Magnetometer  Anomaly  Map  of  Area  1C 
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Figure  14  -  Magnetometer  Anomaly  Map  of  Area  IB 
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Figure  15  -  Magnetometer  Anomaly  Map  of  Area  lA 
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Figure  16  -  Magnetometer  Anomaly  Map  of  Area  BT4  Center 

Three  total  eoverage  areas  were  surveyed  in  the  vieinity  of  Target  3,  loeated  in  the  northern  part 
of  the  WAA  demonstration  site.  The  three  total  coverage  areas  are  labeled  Area  3A,  3B,  and  3C 
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with  Area  3 A  being  the  closest  to  Target  3.  Figure  17,  Figure  18,  and  Figure  19  present  the 
magnetometer  data  anomaly  maps  for  Areas  3A,  3B,  3C  respectively. 
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Figure  17  -  Magnetometer  Anomaly  Map  of  Area  3  A 
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Figure  18  -  Magnetometer  Anomaly  Map  of  Area  3B 
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Figure  19  -  Magnetometer  Anomaly  Map  of  Area  3C 

Two  total  coverage  areas  were  surveyed  in  the  vicinity  of  the  Suspected  75mm  Range  area  of 
interest,  located  in  the  northeastern  portion  of  the  WAA  demonstration  site.  The  two  total 
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coverage  areas  are  labeled  Area  2A  and  2B  with  Area  2A  located  in  the  northwestern  corner  of 
the  Suspected  75mm  Range  area  of  interest  and  Area  2B  located  in  the  southeastern  corner. 
Figure  20  and  Figure  21  present  the  magnetometer  data  anomaly  maps  for  Areas  2A  and  2B 
respeetively.  Area  2A  is  split  vertically  on  the  western  side  by  a  barbed  wire  fence  and  cattle 
guard  on  the  road.  The  survey  was  stopped  several  swath  widths  on  either  side  of  the  fence  to 
limit  the  impact  of  the  fence  on  the  data  collected.  The  southeastern  portion  of  Area  2A  also  had 
a  large  number  of  small  trees  and  eactus  whieh  resulted  in  small  areas  where  data  could  not  be 
collected. 
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Figure  20  -  Magnetometer  Anomaly  Map  of  Area  2A 
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Figure  21  -  Magnetometer  Anomaly  Map  of  Area  2B 
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3,1,6  Anomaly  Density  Falloff  Analysis  for  Targets  3  and  4 

One  intention  of  the  total  coverage  surveys  conducted  in  Areas  1  (Target  4)  and  3  (Target  3)  was 
to  map  the  anomaly  density  falloff  as  a  function  of  distance  from  the  target.  Once  the  total 
coverage  data  had  been  collected  and  analyzed  in  the  MTADS  DAS,  the  data  was  divided  into 
30m  X  30m  cells  in  an  East  /  West  radial  leading  from  the  center  of  each  target  in  an  easterly 
direction.  Figure  22  depicts  the  total  coverage  plan  for  Area  3  (Target  3).  The  red  diamond 
indicates  the  CSM  vO  center  of  the  Target  3  target  circle.  The  red  line  indicates  the  path  of  the 
30m  X  30m  cells.  The  blue  rectangles  represent  the  planned  locations  of  the  total  coverage  areas. 
Some  modification  to  the  area  locations  were  made  in  the  field  and  the  small  mismatch  between 
the  red  line  and  blue  rectangles  reflects  this.  The  number  of  anomalies  in  each  cell  was  counted 
and  is  shown  in  Figure  23.  Assuming  that  the  anomaly  density  around  a  target  falls  off 
according  to  a  normal  distribution,  the  results  can  be  fit  to  a  normal  distribution  with  a  persistent 
background  value.  Such  a  fit  is  shown  in  Figure  23  as  a  solid  line.  If  the  center  of  the 
distribution  is  fixed  at  the  center  of  the  CSM  vO  target  circle,  the  resulting  background  value  is 
1 .8  anomalies  per  cell.  Allowing  the  center  to  float  yields  the  same  background  value  and 
displaces  the  center  position  to  -38m  along  the  radial. 


616500  617000  617500  618000  618500 


Figure  22  -  Total  Coverage  Plan  for  Area  3  (Target  3).  The  planned  total  coverage  survey  areas 
are  shown  in  blue,  the  Target  3  target  circle  from  CSM  vO  is  shown  in  dark  purple  and  the  ASR 
target  outlines  are  shown  in  pink.  The  red  diamond  indicates  the  center  of  the  Target  3  target 
circle.  The  red  line  indicates  the  swath  selected  for  the  radial  analysis. 
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Figure  23  -  Number  of  anomalies  per  30m  x  30m  cell  as  a  function  of  radial  distance  from  the 
CSM  vO  T3  target  circle  center.  The  solid  line  is  the  results  of  a  fit  to  a  normal  distribution  with 
a  persistent  background  value  of  1.8  anomalies  /  cell. 

Figure  24  depicts  the  total  coverage  plan  for  Area  1  (Target  4).  The  red  diamond  indicates  the 
AMT  ADS  magnetometer  survey  center  of  Target  4.  The  red  line  indicates  the  swath  of  the  30m 
X  30m  cells.  The  blue  rectangles  represent  the  planned  locations  of  the  total  coverage  areas. 
Some  modification  to  the  area  locations  were  made  in  the  field  and  the  small  mismatch  between 
the  red  line  and  blue  rectangles  reflects  this.  The  number  of  anomalies  in  each  cell  was  counted 
and  is  shown  in  Figure  25.  Assuming  that  the  anomaly  density  around  a  target  falls  off 
according  to  a  normal  distribution,  the  results  can  be  fit  to  a  normal  distribution  with  a  persistent 
background  value.  Such  a  fit  is  shown  in  Figure  25.  If  the  center  of  the  distribution  is  fixed  at 
the  center  of  the  AMTADS  survey,  the  resulting  background  value  is  1 .4  anomalies  per  cell. 
Table  3-7  tabulates  these  results  along  with  the  overall  anomaly  densities  for  the  Area  2 
(Suspected  75mm  Range)  and  the  Simmons  Area  for  comparison.  The  background  anomaly 
density  1200  -  2000  m  from  the  center  of  the  two  Targets  remains  higher  than  any  other  area 
subjected  to  a  total  coverage  survey. 
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Table  3-7  -  Background  Anomaly  Densities  for  Total  Coverage  Areas  1,  2,  3,  and  the  Simmons 
Area 


Area 

Anomalies  /  30m  x  30m  cell 

Area  1  (Target  4) 

1.4 

Area  2A  (Suspected  75mm  Range) 

1.1 

Area  2B  (Suspected  75mm  Range) 

0.5 

Area  3  (Target  3) 

1.8 

Simmons  Area 
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Figure  24  -  Total  Coverage  Plan  for  Target  4.  The  planned  total  coverage  survey  areas  are 
shown  in  blue,  the  Target  4  target  circle  from  CSM  vO  and  the  ASR  target  outline  are  shown  in 
dark  brown.  The  red  diamond  indicates  the  center  of  Target  3  as  reported  from  the  AMTADS 
magnetometer  data  collected  by  Sky  Research.  The  red  line  indicates  the  swath  selected  for  the 
radial  analysis. 
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Radial  Distance  from  AMTADS  T4  Center  (m) 

Figure  25  -  Number  of  anomalies  per  30m  x  30m  eell  as  a  function  of  radial  distance  from  the 
AMTADS  T4  center.  The  solid  line  is  the  results  of  a  fit  to  a  normal  distribution  with  a 
persistent  background  value  of  1 .4  anomalies  /  cell. 

3.1.7  Section  9  Farm  House  Area  of  Interest 

At  the  request  of  the  ESTCP  Program  Office,  the  transect  survey  data  from  an  identified  area  of 
interest  in  the  central  region  of  the  WAA  demonstration  site  was  loaded  into  the  MTADS  DAS 
software  and  analyzed.  The  permit  holders  were  queried  about  this  area  of  interest  and  reported 
that  it  the  location  of  an  old  homestead  with  a  remaining  structure.  The  tow  vehicle  operator 
reported  a  large  quantity  of  surface  barb-wire  in  the  area  of  the  structure.  The  data  from  Section 
9,  N/S  transects  11-14  and  Area  27,  EAV  transects  1-5,  also  located  in  Section  9  were  analyzed 
and  a  composite  magnetometer  anomaly  map  is  shown  in  Eigure  26.  Eor  the  purposes  of  this 
report  the  area  will  be  referred  to  as  the  Section  9  Earm  House  area  of  interest.  188  anomalies 
were  identified  and  76  of  those  anomalies  were  reasonably  fit  to  the  resident  dipole  model 
(coherence  >  0.85).  The  anomalies  which  could  be  fit  using  the  dipole  model  are  reported  with 
all  parameters  in  the  anomaly  list  on  the  attached  DVD.  The  remaining  anomalies  did  not  fit  well 
enough  to  the  dipole  model  and  only  positions  and  often  a  comment  are  reported. 
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Figure  26  -  Magnetometer  Anomaly  Map  of  Seetion  9  Farm  House  Area  of  Interest 

3,1,8  Calibration  Items 

As  mentioned  in  Section  3.1.1,  a  calibration  strip  of  munitions  and  munitions  stimulants  were 
emplaced  between  the  WAA  Base  Camp  and  the  northern  boundary  of  the  WAA  Demonstration 
site.  Additionally  a  16-lb  shotput  was  emplaced  near  the  Auxiliary  Base  Camp  during  use. 

Table  3-3  gives  a  schedule  of  the  emplaced  items  and  parameters  (i.e.  depth  and  orientation). 
Figure  27  shows  a  magnetometer  anomaly  map  of  the  calibration  strip.  Figure  28  shows  a 
magnetometer  anomaly  map  of  the  calibration  sphere  emplaced  near  the  Auxiliary  Base  Camp. 
For  both  figures,  the  midpoint  positions  of  the  emplaced  items,  as  determined  by  RTK  GPS 
waypointing,  are  shown  as  open  circles. 


36 


616440 


616460 


616480 


Figure  27  -  Magnetometer  anomaly  map  of  the  calibration  strip  emplaced  between  the  WAA 
Base  Camp  and  the  WAA  Demonstration  site  at  Pueblo  PBR#2 
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Figure  28  -  Magnetometer  anomaly  map  of  the  calibration  sphere  emplaced  near  the  Auxiliary 
Base  Camp  located  south  of  the  WAA  demonstration  site  at  Pueblo  PBR#2 

Each  survey  day  commenced  with  collection  of  a  5-6  minute  static  survey  after  the  sensors  had 
been  warmed  up  and  RTK  GPS  was  established.  After  the  static  survey,  the  calibration  items 
were  surveyed,  the  calibration  strip  when  operations  were  being  conducted  out  of  the  WAA  Base 
Camp  and  the  calibration  sphere  when  operations  were  being  conducted  out  of  the  Auxiliary 
Base  Camp.  At  the  end  of  the  survey  day,  the  appropriate  calibration  area  was  surveyed  again. 

In  some  cases  the  survey  day  would  begin  at  one  base  camp  and  end  at  the  other.  In  this 
situation,  the  nearest  calibration  area  was  surveyed.  To  evaluate  the  data  from  the  calibration 
items,  the  signal  to  noise  ratio  (SNR)  was  determined.  The  peak  amplitude  (analytic  signal)  of 
each  calibration  item  emplaced  in  the  calibration  strip  was  determined  for  each  data  set  using  the 
automated  anomaly  picking  routines  to  determine  a  Signal  value  for  each  item.  A  sub-area 
within  the  bounds  of  all  surveys,  identified  to  be  relatively  free  of  anomalies,  was  used  for  each 
data  set  to  extract  a  small  area  of  the  magnetometer  data  and  convert  it  to  analytic  signal.  The 
analytic  signal  standard  deviation  (la)  was  then  calculated  for  the  sub-area  and  that  value  was 
used  as  the  Noise  value  for  each  survey.  The  bounds  of  the  background  areas  for  the  calibration 
strip  and  calibration  sphere  are  given  in  Table  3-8. 

Table  3-8  -  Boundaries  of  Background  Areas  Selected  for  Noise  Calculation 


Boundary 

Easting  (UTM,m) 

Northing  (UTM,  m) 

Calibration  Strip 

Minimum 

618,445 

4,178,771 

Maximum 

618,457 

4,178,785 

Calibration  Sphere 

Minimum 

618,353 

4,168,710 

Maximum 

618,360 

4,168,717 

The  aggregate  values  of  the  SNR  for  all  data  sets  for  the  calibration  strip  items  (average  and 
standard  deviation  (la)  are  tabulated  in  Table  3-9. 
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Survey  File 

Figure  29  and  Figure  30  plot  the  SNR  values  for  the  155mm  Projectile  #2  and  the  60mm  Mortar 
#1  for  all  data  sets  in  a  pseudo-time  series  as  examples.  The  155mm  Projectile  #2  had  the  largest 
SNR  (and  analytic  signal)  values.  The  60mm  Mortar  #1  values  were  approximately  1/10  the 
155mm  values  and  represent  the  smaller  values  measured.  The  solid  line  indicates  the  aggregate 
average  and  the  dashed  lines  indicate  ala  envelope. 

Table  3-9  -  Signal  to  Noise  (SNR)  Aggregate  Values  for  Calibration  Strip  Items 


ID 

Easting 
(UTM,  m) 

Northing 

(UTM,m) 

Depth  (cm) 

Orientation 
(deg.,  from 
North) 

Avg.  SNR 

SNR  SD  (la) 

Sphere  #1  (Driver  Side) 

616434.500 

4178732.403 

0 

N/A 

623.09 

135.43 

Sphere  #2  (Passenger  Side) 

616435.459 

4178732.109 

0 

N/A 

356.24 

48.86 

155mm  Projectile  #2 

616441.180 

4178749.703 

35 

35 

1914.52 

244.19 

60  mm  Mortar  #2 

616447.267 

4178768.742 

30 

46 

18.19 

2.35 

105mm  Projectile  #2 

616453.828 

4178787.621 

60 

44 

66.26 

7.62 

105mm  Projectile  #1 

616459.639 

4178806.967 

45 

178 

231.75 

31.08 

81mm  Mortar  #2 

616465.341 

4178825.771 

43 

69 

53.03 

6.47 

81mm  Mortar  #1 

616469.792 

4178839.941 

25 

20 

90.97 

9.72 

155mm  Projectile  #1 

616474.350 

4178854.198 

50 

46 

536.95 

61.99 

60  mm  Mortar  #1 

616478.260 

4178868.186 

10 

148 

107.68 

14.20 

37mm  Sim  #2 

616481.220 

4178877.825 

10 

57 

139.69 

18.25 

37mm  Sim  #1 

616484.096 

4178887.419 

5 

160 

255.47 

37.02 
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Figure  29  -  SNR  data  runs  for  155mm  Projectile  #2.  The  result  for  each  data  set  is  shown  in 
order  of  acquisition.  The  horizontal  axis  is  survey  file  number.  The  solid  line  represents  the 
aggregate  average  SNR  and  the  dashed  lines  represent  ala  envelope. 
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Figure  30  -  SNR  data  runs  for  60mm  Projectile  #1 .  The  result  for  each  data  set  is  shown  in  order 
of  acquisition.  The  horizontal  axis  is  survey  file  number.  The  solid  line  represents  the  aggregate 
average  SNR  and  the  dashed  lines  represent  ala  envelope. 
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The  aggregate  values  of  the  analytie  signal  alone  for  all  data  sets  for  the  calibration  strip  items 
(average  and  standard  deviation  (la)  are  tabulated  in  Table  3-10.  Figure  31  and  Figure  32  plot 
the  Signal  values  for  the  155mm  Projectile  #2  and  the  60mm  Mortar  #1  for  all  data  sets  in  a 
pseudo-time  series  as  examples.  The  155mm  Projectile  #2  had  the  largest  signal  values.  The 
60mm  Mortar  #1  values  were  approximately  1/10  the  155mm  values  and  represent  the  smaller 
values  measured.  The  solid  line  indicates  the  aggregate  average  and  the  dashed  lines  indicate  a 
1g  envelope.  Comparing  the  SNR  and  analytic  signal  values  indicates  that  the  analytic  signal 
values  are  approximately  twice  that  of  the  SNR  (typical  Pueblo  PBR#2  Noise  values  being  on 
the  order  of  2  nT/m)  and  have  a  standard  deviation  that  is  one  half  the  SNR  value.  This  indicates 
that  measuring  the  peak  amplitude  of  the  item  analytic  signal  is  not  the  limiting  factor  in  the 
SNR  calculation  in  this  case. 

Table  3-10  -  Signal  Aggregate  Values  for  Calibration  Strip  Items 


ID 

Easting 
(UTM,  m) 

Northing 

(UTM,m) 

Depth  (cm) 

Orientation 
(deg.,  from 
North) 

Avg.  Signal 
(nT/m) 

Signal  SD 
(la,  nT/m) 

Sphere  #1  (Driver  Side) 

616434.500 

4178732.403 

0 

N/A 

1064.56 

203.42 

Sphere  #2  (Passenger  Side) 

616435.459 

4178732.109 

0 

N/A 

645.75 

47.30 

155mm  Projeetile  #2 

616441.180 

4178749.703 

35 

35 

3461.28 

153.18 

60  mm  Mortar  #2 

616447.267 

4178768.742 

30 

46 

32.90 

1.96 

105mm  Projeetile  #2 

616453.828 

4178787.621 

60 

44 

119.92 

4.38 

105mm  Projeetile  #1 

616459.639 

4178806.967 

45 

178 

418.80 

21.70 

81mm  Mortar  #2 

616465.341 

4178825.771 

43 

69 

95.98 

5.02 

81mm  Mortar  #1 

616469.792 

4178839.941 

25 

20 

165.35 

14.69 

155mm  Projeetile  #1 

616474.350 

4178854.198 

50 

46 

971.91 

37.72 

60  mm  Mortar  #1 

616478.260 

4178868.186 

10 

148 

195.48 

20.91 

37mm  Sim  #2 

616481.220 

4178877.825 

10 

57 

253.46 

24.75 

37mm  Sim  #1 

616484.096 

4178887.419 

5 

160 

463.54 

53.27 
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Figure  31  -  Analytic  signal  data  runs  for  155mm  Projectile  #2.  The  result  for  each  data  set  is 
shown  in  order  of  acquisition.  The  horizontal  axis  is  survey  file  number.  The  solid  line 
represents  the  aggregate  average  analytic  signal  and  the  dashed  lines  represent  ala  envelope. 
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Figure  32  -  Analytic  signal  data  runs  for  60mm  Mortar  #1 .  The  result  for  each  data  set  is  shown 
in  order  of  acquisition.  The  horizontal  axis  is  survey  file  number.  The  solid  line  represents  the 
aggregate  average  analytic  signal  and  the  dashed  lines  represent  ala  envelope. 
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A  similar  analysis  was  conducted  for  the  metal  sphere  emplaced  near  the  Auxiliary  Base  Camp. 
The  results  and  the  aggregate  values  of  the  SNR  and  the  analytic  signal  for  all  data  sets  for  the 
calibration  sphere  (average  and  standard  deviation  (la)  are  shown  in  Figure  33  and  Figure  34. 
The  aggregate  values  for  the  metal  sphere  were  an  SNR  of  191.9  ±  79.1  and  a  peak  analytic 
signal  amplitude  of  857.7  ±  235.9  nT/m.  If  the  analysis  is  restricted  to  the  morning  calibration 
runs  only,  the  standard  deviation  reduces  by  approximately  one-half  to  values  of  185.6  ±  44.5  for 
the  SNR  and  919.6  ±  121.5  nT/m  for  the  analytie  signal.  The  tow  vehicle  batteries  were  charged 
by  a  5  kW  rather  than  a  65  kW-range  generator  at  the  Auxiliary  Base  Camp.  This  may  show  a 
limit  of  the  system’s  daily  operational  endurance  when  charging  the  batteries  overnight  using  the 
smaller  generator  sinee  this  phenomenon  was  not  observed  in  the  data  collect  when  operating  out 
of  the  WAA  Base  Camp. 
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Figure  33  -  SNR  data  runs  for  the  metal  sphere  emplaced  near  the  Auxiliary  Base  Camp.  The 
result  for  each  data  set  is  shown  in  order  of  aequisition.  The  horizontal  axis  is  survey  file 
number.  The  solid  line  represents  the  aggregate  average  SNR  and  the  dashed  lines  represent  a 
la  envelope. 
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Figure  34  -  Analytic  signal  data  runs  for  the  metal  sphere  emplaced  near  the  Auxiliary  Base 
Camp.  The  result  for  each  data  set  is  shown  in  order  of  acquisition.  The  horizontal  axis  is 
survey  file  number.  The  solid  line  represents  the  aggregate  average  analytic  signal  and  the 
dashed  lines  represent  ala  envelope. 

As  stated  above,  static  tests  of  the  sensor  platform  were  conducted  at  the  beginning  of  each 
survey  day.  After  a  period  for  system  warm  up  (approximately  15  minutes),  walk-around 
preventative  maintenance  inspections,  and  the  establishment  of  RTK  GPS,  the  tow  vehicle  was 
driven  to  the  appropriate  calibration  area.  A  data  set  was  collected  for  at  least  5-6  minutes  while 
the  vehicle  was  kept  stationary.  The  2-D  positioning  variation  was  evaluated  by  computing  the 
standard  deviation  of  both  the  northing  and  easting  components  of  the  position  data  for  the  entire 
period  and  combining  them  as  the  square  root  of  the  sum  of  the  squares.  The  standard  deviation 
for  the  demedianed  magnetometer  data  from  each  sensor  was  computed  and  the  arithmetic  mean 
computed  was  for  the  data  set.  In  occasional  cases,  an  obvious  artifact  was  present  in  the  data 
(e.g.  a  vehicle  pulls  up  along  side  the  tow  vehicle  unannounced)  and  distorts  a  portion  of  the 
static  run.  In  these  cases,  only  the  unperturbed  data  was  used  with  a  minimum  of  3  minutes  total 
data.  The  standard  deviation  of  both  the  positioning  and  sensor  data  for  all  data  sets  was 
computed.  The  results  are  shown  in  the  follow  pseudo-time  series  figures.  Figure  35  and  Figure 
36  show  the  positioning  and  magnetometer  variations  for  the  static  tests  conducted  at  the 
southern  end  of  the  calibration  strip.  Figure  37  and  Figure  38  display  the  results  of  a  similar 
analysis  for  the  data  sets  collected  at  the  western  edge  of  the  survey  path  used  for  the  calibration 
sphere  located  near  the  Auxiliary  Base  Camp. 

The  magnetometer  values  shown  in  Figure  38  indicate  a  significant  increase  in  the  overall 
average  magnetometer  sensor  noise  level  for  Julian  dates  05277  and  later.  The  first  three-week 
portion  of  the  demonstration  depleted  the  reserve  of  spare  magnetometers.  The  manufacturer 
had  promised  to  deliver  the  repaired  spares  prior  to  the  commencement  of  the  second  portion  of 
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the  survey.  The  return  of  the  sensors  was  delayed  due  to  difficulties  in  tuning  the  sensors  at  the 
factory.  The  repaired  sensors  arrived  after  the  survey  was  stopped  on  October  7  for  unscheduled 
repairs  to  the  vehicle  tow  vehicle  engine.  A  consequence  was  that  the  surveys  were  conducted 
with  one  sensor  operating  at  a  noise  level  approximately  10  times  the  norm.  The  data  from  this 
sensor  was  only  used  for  further  analysis  when  the  SNR  was  10  or  greater.  For  the  third  portion 
of  the  survey  the  spare  sensors  were  deployed  as  required.  An  intermittent  problem  with  the 
counter  board  for  one  outboard  sensor  developed  during  the  third  portion  of  the  survey.  A  part 
was  required  to  troubleshoot  /  repair  this  issue  that  is  not  typically  kept  in  the  spares  collection. 
The  data  from  this  sensor  was  not  used  for  further  analyses  for  data  sets  when  the  problem 
manifest  itself  while  the  part  was  being  delivered.  The  part  arrived  on  October  21  and  the 
problem  did  not  manifest  itself  again.  The  combination  of  these  issues  results  in  the  increase 
sensor  measurement  variation  indicated  in  Figure  38.  In  the  cases  where  all  sensors  appeared  to 
be  working  within  limits  and  yet  the  static  variation  is  higher  than  expected,  a  long  wavelength 
variation  in  the  sensors  could  be  observed  in  the  data  but  not  by  the  operator  in  the  field.  It  is 
possible  that  the  source  of  this  issue  was  insufficient  warm  up  time  in  the  morning  but  the  issue 
is  still  under  investigation  and  only  occurred  during  operations  out  of  the  Auxiliary  Base  Camp 
with  the  smaller  generator.  Table  3-11  summarizes  the  static  test  data  results. 

Table  3-1 1  -  Static  Test  Data  Results 


Calibration  Area 

Result  Type 

Value 

North 

2-D  Position 

0.42  ±  0.14  cm 

Magnetometer 

0.89  ±  0.97  nT 

South 

2-D  Position 

0.44  ±0.10  cm 

Magnetometer 

0.67  ±  0.72  nT 
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Figure  35  -  Positional  variation  data  runs  for  static  data  collected  at  the  calibration  strip.  The 
horizontal  axis  is  survey  file  number.  The  solid  line  represents  the  aggregate  average  positional 
variation  and  the  dashed  lines  represent  ala  envelope. 


Survey  File 


Survey  File 


Figure  36  -  Overall  magnetometer  (all  sensors)  variation  data  runs  for  static  data  collected  at  the 
calibration  strip.  The  horizontal  axis  is  survey  file  number.  The  solid  line  represents  the 
aggregate  average  sensor  variation  and  the  dashed  lines  represent  ala  envelope. 
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Figure  37  -  Positional  variation  data  runs  for  static  data  collected  at  the  calibration  sphere.  The 
horizontal  axis  is  survey  file  number.  The  solid  line  represents  the  aggregate  average  positional 
variation  and  the  dashed  lines  represent  ala  envelope. 
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Figure  38  -  Overall  magnetometer  (all  sensors)  variation  data  runs  for  static  data  collected  at  the 
calibration  sphere.  The  horizontal  axis  is  survey  file  number.  The  solid  line  represents  the 
aggregate  average  sensor  variation  and  the  dashed  lines  represent  ala  envelope. 
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3.1,9  Demobilization 


At  the  end  of  field  operations,  all  equipment,  materials,  and  supplies  were  repaeked  on  the  53’ 
trailer  and  seeured.  Harris  Transportation  Company,  a  government  eontraet  transportation  firm 
transported  the  trailer  from  the  site  to  the  MTADS  home  base  at  ART  Blossom  Point,  Weleome, 
MD.  Onee  all  personnel  had  departed  the  site  the  loeal  vendors  were  eontaeted  to  arrange  for  the 
generator  and  fuel  storage  to  be  removed.  In  eooperation  with  the  land  owner,  Tead  Russell,  the 
office  and  garage  trailers  and  the  portable  toilet  at  the  WAA  Base  Camp  are  being  left  in  place 
for  the  validation  effort  currently  being  planned  for  Spring,  2006  by  the  WAA  Pilot  Project.  The 
second  garage  trailer,  located  at  the  south  end  of  the  site,  was  also  removed  from  its  location  on 
U.S.  Forest  Service  land. 
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Appendix  A.  Analytical  Methods  to  Support  the  Experimental  Design 


To  facilitate  the  scope  and  tempo  of  the  WAA  pilot  project,  the  typical  MTADS  man-in-the-loop 
analysis  of  total  field  magnetometer  data  was  not  going  to  fit  the  bill.  The  proeess  is  too  time 
eonsuming  and  provides  far  rieher  results  than  is  initially  required  by  the  planning  elements  of 
the  project  on  a  daily  time  seale.  An  alternate  approaeh  was  required.  Briefly,  the  loeated  total 
magnetie  field  data  (nT)  are  eolleeted  as  normal  for  an  MTADS  survey.  The  total  field  data  are 
eonverted  to  analytie  signal  (AS,  nT/m)  where  the  analytie  signal  is  ealeulated  from  the  squares 
of  the  derivatives  in  the  x,  y,  and  z  direetions: 


The  z  derivatives  ean  be  estimated  from  the  total  field  2-D  data  using  either  FFT  or  eonvolution 
methods.  For  this  work,  the  FFT  method  was  used.  All  proeessing  work  was  done  in  the 
eommereial  geophysics  software  paekage.  Oasis  montaj  and  requires  the  advaneed  1-D  filters 
available  as  part  of  the  Geophysies  add-on  paekage  (and  others).  A  data  grid  is  ealeulated  for  the 
total  field  data  and  this  involves  interpolating  the  aetual  data  to  a  regularly  spaeed  grid  of  a  given 
eell  size.  The  algorithms  whieh  eompute  the  AS  require  a  regularly  spaeed  input  data  grid  to 
work  efficiently  and  output  an  AS  data  grid.  The  GRIDPEAK  GX  is  then  used  to  extract  all 
peaks  in  the  AS  grid  whieh  are  above  a  given  threshold.  This  Appendix  deseribes  initial 
evaluation  work  done  on  existing  data  sets  to  validate  this  method  and  provide  initial  estimates 
of  the  optimized  parameters  required. 

Total  field  magnetometer  and  AS  grids  were  ealeulated  from  an  existing  data  set  eolleeted  using 
the  MTADS  platform  at  our  home  base  in  Blossom  Point,  MD  at  0.125,  0.25,  0.50,  and  1.00  m 
grid  eell  sizes.  The  GRIDPEAK. GX  was  then  run  using  a  set  of  Cut-Off  threshold  values 
(required  peak  amplitude  to  be  eounted  as  a  peak)  and  a  fixed  number  of  smoothing  filter  passes 
(six)  on  eaeh  AS  grid.  Six  passes  was  found  to  be  the  minimum  to  get  reasonable  results.  Using 
fewer  passes  result  in  multiple  pieks  for  the  same  actual  peak  in  the  grid.  The  following  four 
graphs  (Eigures  A1  -  A4)  show  the  #  of  peaks  pieked  for  a  given  Cut-Off  value  for  eaeh  AS  grid 
eell  size.  The  red  line  indieates  the  number  of  emplaeed  targets  (all  types)  in  the  BP  test  field 
(61). 

Eaeh  curve  shows  a  'knee',  below  whieh  one  picks  up  non-emplaced  targets  faster  than  emplaeed 
targets.  The  0.25m  studies  eaptured  51  to  57  of  the  61  emplaeed  targets  with  a  range  of  non- 
emplaeed  targets  pieked.  A  majority  of  the  undeteeted  emplaeed  targets  were  in  the  Clutter 
eategory  and  the  exelusion  of  these  objects  is  not  necessarily  a  failing  of  the  method. 
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Figure  A1  -  Blossom  Point  GRIDPEAK  Results  for  a  0.125m  grid  eell  size 


Figure  A2  -  Blossom  Point  GRIDPEAK  Results  for  a  0.25m  grid  eell  size 
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Figure  A3  -  Blossom  Point  GRIDPEAK  Results  for  a  0.50m  grid  eell  size 


Figure  A4  -  Blossom  Point  GRIDPEAK  Results  for  a  1.00m  grid  eell  size 


Eor  eight  points  from  the  above  studies  (predominately  from  the  0.25m  grid  eell  size  results),  the 
pieked  anomalies  were  matched  to  an  emplaced  target.  If  the  pick  was  not  clearly  associated 
with  an  emplaced  target,  it  was  marked  as  a  Ealse  Alarm  (EA)  since  the  Blossom  Point  Test  Eield 
is  relatively  clean  after  repeated  clearances  and  should  have  no  non-emplaced  targets.  Based  on 
these  results,  values  for  probability  of  detection  (Pd),  false  alarm  rate  (PAR,  PA's  /  hectare), 
average  miss  distance,  and  the  standard  deviation  of  the  miss  distance  were  calculated. 
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Table  A-1  gives  the  results.  A  ROC  eurve  was  eonstructed  and  is  shown  as  Figure  A5  for  all 
analyzed  cases  with  the  0.25m  cases  plotted  in  red. 

Table  A-1  -  Results  for  emplaced  targets  at  Blossom  Point  for  various  parameters 


Grid  Cell  Size  (m) 

Cutoff  Threshold  (nT\m) 

#  of  Targets  Picked 

#  of  Emplaced  Targets  Picked 

Pd 

#  ot  1  argets  Picked  Not 
Emplaced  Targets 

FAR  («  FA/ 
Hectare) 

Average  Miss 
Distance  (m) 

Std.  Dev 
(1s)  (m) 

0.125 

25 

120 

57 

0.934 

63 

223 

0.119 

0.079 

0.125 

50 

80 

57 

0.934 

23 

82 

0.119 

0.079 

0.250 

25 

70 

57 

0.934 

13 

46 

0.132 

0.080 

0.250 

37 

66 

56 

0.918 

10 

35 

0.133 

0.081 

0.250 

50 

62 

55 

0.902 

7 

25 

0.131 

0.080 

0.250 

75 

52 

51 

0.836 

1 

4 

0.136 

0.082 

0.500 

25 

11 

11 

0.180 

0 

0 

0.419 

0.396 

0.500 

50 

5 

5 

0.080 

0 

0 

0.370 

0.289 

Figure  A5  -  ROC  curve  for  emplaced  target  comparisons 

As  stated  above,  the  BP  Test  Field  is  a  relatively  small  site  with  a  limited  number  of  emplaced 
targets  and  is  relatively  clean.  To  further  test  the  developed  method,  existing  MTADS 
magnetometer  data  collected  at  the  Standardized  UXO  Technology  Demonstration  Sites  located 
at  the  Aberdeen  and  Yuma  Proving  Grounds  was  considered.  Using  0.25m  grid  spacing  and  the 
same  range  of  Cut-Off  threshold  values,  the  #  of  picks  vs.  Cut-Off  Threshold  value  for  these  two 
sites  was  determined.  Plots  are  included  showing  ATC  (Figure  A6)  and  YTC  (Figure  A7) 
separately  and  then  together  with  the  BP  results  for  0.25  m  (Figure  A8)  are  shown.  Notice  the 
strong  similarity  between  the  ATC  and  YTC  results.  The  similarity  to  the  BP  is  also  good 
considering  how  much  smaller/cleaner  the  BP  field  is.  These  results  support  the  general 
applicability  of  the  BP  analysis  to  ATC  and  YTC  (if  no  further).  No  count  of  the  actual  number 
of  emplaced  targets  is  available  for  these  sites,  so  no  reference  lines  are  plotted. 
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Figure  A6  -  ATC  GRIDPEAK  Results  for  a  0.25m  grid  cell  size 


Figure  A7  -  YTC  GRIDPEAK  Results  for  a  0.25m  grid  cell  size 
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Figure  A8  -  YTC,  ATC,  and  BP  GRIDPEAK  Results  for  a  0.25m  grid  eell  size 

The  0.25m  grid  eell  size  appears  to  be  a  promising  combination  of  sensitivity  (from  this  work) 
and  computational  speed.  Timed  trials  conducted  on  a  large  data  set  found  that  the  typical  grid 
cell  sized  used  in  total  coverage  surveys  (0.125  m)  is  computationally  impractical  for  transect 
style  data  with  long  (~6  km)  down  track  distances  and  narrow  cross  track  width  (2  km).  If  the 
computations  were  completed,  the  time  required  was  on  the  order  of  4  hours  per  one  hour  of 
survey. 

In  conclusion,  the  following  three  points  can  be  drawn: 

1)  The  0.25m  grid  cell  size  appears  to  be  a  promising  combination  of  sensitivity  and 
computational  speed. 

2)  A  generally  applicable  guideline  of  25  to  100  nT/m  for  the  AS  cut-off  threshold  value  is 
beginning  to  appear  and  holds  on  three  sites. 

3)  The  choice  of  AS  cut-off  threshold  value  will  be  driven  by  the  relative  tolerance  for  detection 
of  smaller  clutter  items  versus  desired  sensitivity. 

Ongoing  efforts  which  will  feed  directly  into  this  project  include  the  analysis  of  other  data  sets 
and  further  exploration  of  the  available  parameter  space.  One  other  data  set  from  the  MTADS 
BBR  2002  vehicular  survey  exists.  This  data  set  represents  a  plains-like  geology  which  may  be 
relevant  to  the  PBR#2  site.  The  large  positional  errors  reported  for  the  0.5  and  1.0  m  grid  cell 
sizes  leave  room  for  improvement.  Alternate  smoothing  filter  settings  may  prove  more  effective 
for  the  larger  grid  cell  sizes.  This  will  be  an  area  of  continued  effort. 
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Appendix  B.  Data  Storage  and  Archiving  Procedures 

B,1  Data  Formats 

Each  survey  file  set  eontains  10  files  which  constitute  the  ‘raw  data’.  The  file  name  strueture  is 
YYDDDFFF.DevieeType.DevieeAlias;  where  YY  is  the  2-digit  year,  DDD  is  the  "Julian"  day, 
or  day  in  the  year,  and  FFF  is  the  flight  number  starting  with  001 .  In  the  following  example,  the 
data  was  taken  on  the  210th  day  of  2002,  flight  number  4. 

02210004.Survey.822A.822A_l 

02210004.Survey.822A.  822A_2 

02210004.  Survey .  GPS  .NMEA 

022 10004.  Survey. SerialDevice.UTC 

022 10004. Survey. PpsDeviee.PPS 

022 10004. Survey. TriggerDeviee. Trigger 

022 10004. Survey. FineNumber 

02210004.Survey 

022 10004.  Survey  .page 

022 10004.  Survey  .loginfo  1  .txt 

Eaeh  data  line  is  time  stamped  with  the  PC  system  eloek  to  allow  syncronization  between  files 

YYDDDFFF. Survey. LineNumber  -  start  and  stop  time  of  eaeh  line  in  survey,  typieally  only  one  line  /  file 
YYDDDFFF.Survey.822A.822A_l  -  Output  from  Counter  1  (4  magnetometers),  in  nT  x  lO^'S,  50  Flz. 
YYDDDFFF. Survey. 822 A. 822A_2  -  Output  from  Counter  2  (4  magnetometers),  in  nT  x  10^5,  50  Hz. 
YYDDDFFF. Survey.PpsDeviee.PPS  -  pulse  per  seeond  (PPS)  from  GPS  reeeiver,  1  Hz. 

YYDDDFFF. Survey.GPS.NMEA  -  GPS  output,  Trimble  PTNL,GGK  sentenee  at  5  Hz  (position). 

YYDDDFFF. Survey.TriggerDeviee. Trigger  -  trigger  pulse  to  magnetometers,  50  Hz. 

YYDDDFFF. Survey.SerialDeviee.UTC  -  UTC  time  tag  from  GPS  reeeiver,  "The  time  will  be"  message  for  next 
PPS,  1  Hz. 

The  .Survey,  .Survey. page,  and  .Survey. loginfo*. txt  files  are  setup  information  reeorded  by  the 
data  eolleetion  program  and  contain  no  data  of  use  to  the  user. 


.Survey. FineNumber  files: 

START  LINE  0  12/21/04  12:45:39.523 

STOP  LINE  0  12/21/04  12:59:21.072 


Magnetometer  (.822A)  files: 


dl5289543808d25289567673d35289555967d45289802122 

dl5289545560d25289568728d35289557064d45289803821 

dl5289547878d25289569235d35289557743d45289805162 

dl5289547468d25289568538d35289557255d45289804417 

dl5289546204d25289567936d35289556456d45289802950 

dl5289545018d25289566714d35289556217d45289801466 


10/10/02 

10/10/02 

10/10/02 

10/10/02 

10/10/02 

10/10/02 


14:17:00.508 

14:17:00.528 

14:17:00.548 

14:17:00.568 

14:17:00.588 

14:17:00.608 


First  line: 

dl  -  Sensor  1  ok  -  two  charaeters  of  status  code  /  marker  -  other  two  charaeter  codes  are  possible 
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to  indicate  error  eonditions 

5289543808  -  52895.43808  gamma  or  nT 

d2  -  Sensor  2  ok 

5289567673  -  52895.67673  nT 

d3  -  Sensor  2  ok 

5289555967  -  52895.55967  nT 

d4  -  Sensor  2  ok 

5289802122  -  52898.02122  nT 

10/10/02  -  eomputer  date  stamp  for  reeeipt  of  string  at  eomputer. 

14:17:00.508  -  eomputer  time  stamp  for  reeeipt  of  string  at  eomputer. 

•Survey.PpsDeviee.PPS  fdes: 

PPS  12/21/04  12:45:40.433 

PPS  12/21/04  12:45:41.433 

PPS  12/21/04  12:45:42.433 

. Survey. GPS.NMEA  files: 

$PTNL,GGK, 175017. 00, 122104, 3825. 0 6336634, N, 07706. 2665 6042, W, 3, 07, 2 . 8,EHT- 
25.694,M*7C  12/21/04  12:45:39.470 

Table  C-1  -  PTNL,GGK  Message  Fields^ 


Field 

Meaning 

1 

UTC  of  position  fix 

2 

Date 

3 

Latitude 

4 

Direction  of  Latitude  (N  =  North,  S  =  South) 

5 

Longitude 

6 

Direction  of  Longitude  (E  =  East,  W  =  West) 

7 

GPS  Fix  Quality  (0  =  Invalid,  1,2,3, 4) 

8 

Number  of  Satellites  in  fix 

9 

DOP  of  fix 

10 

Ellipsoidal  height  of  fix 

11 

M:  ellipsoidal  height  is  measured  in  meters 

“  For  further  information,  refer  to  the  Trimble  MS  Series  Operation 
Manual 


. Survey. SerialDeviee.UTC  files: 

UTC  04.12.21  17:50:18  57  12/21/04  12:45:39.645 

UTC  04.12.21  17:50:19  57  12/21/04  12:45:40.646 


Loeated  data  archives  are  ASCII  files  of  the  format: 
For  loeated,  demedianed  magnetometer  data: 
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X  (UTM  Zone  X,  NAD83,  m)  Easting 

Y  (UTM  Zone  X,  NAD83,  m)  Northing 

Z  Height  Above  Ellipsoid  (HAE,  WGS84,  m) 

S  Signal  in  nT 

where  X  is  the  appropriate  UTM  zone  (13N  for  PBR#2,CO) 
for  the  analytie  signal  data,  the  Signal  is  reported  in  nT/m. 

Anomaly  Report  (.Anomaly)  Files: 

Anomaly  Reports  from  Transeet  data  will  be  ASCII  files  of  the  format: 

ID  Fiducial  ID  of  the  anomaly 
X  (UTM  Zone  X,  NAD83,  m)  Easting 

Y  (UTM  Zone  X,  NAD83,  m)  Northing 
S  Analytic  Signal  in  nT\m 

where  X  is  the  appropriate  UTM  zone  (13N  for  PBR#2,CO) 

Course  over  Ground  (.COG)  fdes: 

Corresponding  Course-Over-Ground  (COG)  Reports  for  Transect  data  will  be  ASCII  fdes  of  the 
format: 

X  (UTM  Zone  X,  NAD83,  m)  Easting 

Y  (UTM  Zone  X,  NAD83,  m)  Northing 

GPSTime  UTC  Time  in  seconds  past  midnight 

where  X  is  the  appropriate  UTM  zone  (I3N  for  PBR#2,CO) 

Static  Survey  Archive  (_static.xyz)  files: 

Daily  static  calibration  run  data  will  be  archived  as  geosoft  .XYZ  files  of  the  format: 


X 

(UTM 

Zone  X,  NAD83,  m)  Easting  for 

■  GPS 

antenna 

Y 

(UTM 

Zone  X,  NAD83,  m)  Northing  for  GPS  antenna 

HAE 

(WGS84,  m)  Height  above  Ellipsoid 

for 

GPS  antenna 

Magi 

(nT) 

Demedianed 

magnetometer 

data 

for 

sensor 

1 

Mag2 

(nT) 

Demedianed 

magnetometer 

data 

for 

sensor 

2 

Mag  3 

(nT) 

Demedianed 

magnetometer 

data 

for 

sensor 

3 

Mag  4 

(nT) 

Demedianed 

magnetometer 

data 

for 

sensor 

4 

Mag  5 

(nT) 

Demedianed 

magnetometer 

data 

for 

sensor 

5 

Mag  6 

(nT) 

Demedianed 

magnetometer 

data 

for 

sensor 

6 

Mag  7 

(nT) 

Demedianed 

magnetometer 

data 

for 

sensor 

7 

Mag  8 

(nT) 

Demedianed 

magnetometer 

data 

for 

sensor 

8 

where  X  is  the  appropriate  UTM  zone  (13N  for  PBR#2,CO) 


MTADS  DAS  Target  List  Example 

The  example  is  given  in  ASCII  text  file  format.  Actual  delivery  will  be  in  Excel  Spreadsheet 
format. 

MTADS  TARGET  REPORT 

################### 

Mon  Oct  31  14:00:47  2005 
PROJECT :  PBR2 
SITE:  Area_2A 
SENSOR:  mag 
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SURVEY:  Survey 

PRIMARY  COORDINATES:  UTM=13,  nad83 
################### 

################### 

################### 

################### 

ID,UTM  X  (m)  ,  UTM  Y  (m)  ,  Depth  (in),Size  (m).  Moment  (Amps- 
m2 ), Inclination, Azimuth, Goodness  of  Fit, Comments 
1, 617 608. 50, 417 6876. 99, 0.331, 0.02 8, 0.012 1,26. 70, 30. 32, 0.97 14, 
2, 6177 93. 59, 417 6877. 94, 0.931, 0.1 68, 2. 57 17, 57. 71, 10. 09, 0.9362, 
3,  617799. 14, 417 68 67. 65,  0.844, 0.12  5,  1.04  76,  55. 00, 24. 48, 0.9 964, 
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